Attracting shallow donors: Hydrogen passivation in (AI,Ga,ln)-doped ZnO 
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The hydrogen interstitial and the substitutional Alz n , Gaz n and Inz n are all shallow donors in ZnO and lead 
to n-type conductivity. Although shallow donors are expected to repel each other, we show by first principles 
calculations that in ZnO these shallow donor impurities attract and form a complex, leading to a donor level 
deep in the band gap. This puts a limit on the n-type conductivity of (Al,Ga,In)-doped ZnO in the presence 
of hydrogen. 



The incorporation of hydrogen in a material can 
strongly affect its electronic properties. Furthermore, its 
behavior depends on the host into which it is introduced. 
Mostly, interstitial hydrogen is amphoteric, i.e. it can act 
either as a donor (Hi + ) or an acceptor (Hi~)i. Which of 
these two prevails depends on the Fermi level: it acts as 
a donor in p-type materials and as an acceptor in n-type 
materials. Consequently, it counteracts the conductivity 
caused by extrinsic dopants. It is important to realize 
that this passivation of extrinsic dopants by hydrogen is 
not caused by the formation of a neutral complex. The 
formation of such a neutral complex may just be the con- 
sequence after compensation, i.e. after the trapping of a 
free hole or electron by hydrogen. Afterwards, the posi- 
tively charged H and negatively charged acceptor, or neg- 
atively charged H and positively charged donor attract 
and may form a neutral complex. 

In ZnO however, the behavior of hydrogen is not am- 
photeric 1 . Nominally undoped ZnO shows n-type con- 
ductivity and it was first predicted by first-principles cal- 
culations based on density functional theory^, and later 
experimentally verified 3 -, that interstitial hydrogen acts 
as a shallow donor in ZnO and can cause this n-type 
behavior. 

ZnO is a wide band gap (~3.4 eV at room tempera- 
ture 4 ) semiconductor with many possibilities for techno- 
logical applications as a transparent conductor in solar 
cells 5 , flat panel displays 6 , etc. It is well-known that the 
electronic conductivity of ZnO can be dramatically im- 
proved by doping it with group III elements Al, Ga or 
In 7 -. This makes ZnO an alternative for indium tin ox- 
ide, which shows highest conductivity among transpar- 
ent conducting oxides but is less abundant and therefore 
more expensive. It is then also important to know what 
limits the conductivity of ZnO doped by Al, Ga or In. 
Properties such as formation energy, energy band struc- 
ture and equilibrium geometry of Al, Ga and In doped 
ZnO were studied systematically and it was shown that 
all three dopants are indeed shallow donors and that they 
possess the intrinsic qualities to be good n-type transpar- 
ent conducting oxides 8 . 
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As the concentration becomes high, the interaction be- 
tween defects becomes unavoidable. Because it is very 
difficult to remove hydrogen from the crystal growth envi- 
ronment, it is then also natural to ask what the influence 
is of interstitial hydrogen on Al, Ga or In doped ZnO. 
Recently a number of experiments were reported where 
both H and group III elements are co-doped into ZnO. It 
was shown that the conductivity of ZnO co-doped with 
H and group III elements is better than that doped with 
group III elements up to certain amount of H inclusion 9 . 

Apart from single hydrogen dopants, the influence of 
hydrogen on doped ZnO structures has been studied be- 
fore by first principles calculations. It was found that 
hydrogen passivates nitrogen dopants 10 , silver dopants 11 
and Zn vacancies^ in ZnO, and forms complexes. This is 
not very surprising, as in all these cases, hydrogen, which 
is a (positive) shallow donor, is attracted to a (negative) 
acceptor impurity. This is not unique to hydrogen. Also 
donor-like Asz n and acceptor-like Vz n attract each other 
and form a complex in As doped ZnO 13 . Recently, also 
the attractive interaction between two donor impurities 
was predicted by first principles studies in ZnO, namely 
the deep donor Vo (the oxygen vacancy) and the shal- 
low donor Zn^ (the Zn interstitial). In this Letter, we 
consider the two shallow donors interstitial hydrogen 
and substitutional Xz n with X = Al, Ga or In. 

We performed density functional theory calculations as 
implemented in the Vienna ab initio simulation package 
(VASP) cod o 15 i 16 . In order to examine interactions be- 
tween defects, a large supercell consisting of 192 atoms 
was considered 17 . A single Al, Ga or In atom was chosen 
to substitute a Zn atom, while for the H atom, an exten- 
sive number of positions around the group III element 
were considered. For all these configurations the forma- 
tion energies for different charge states were calculated 
using the method of Ref. 18 . This method was recently 
confirmed to give reliable results for the transition ener- 
gies between different charge states of the impurity, both 
for ZnO with interstitial Hi£ as for ZnO doped with a 
group III element 8 , by comparing the results with hybrid 
functional calculations. 

For all charge states, we find that the lowest formation 
energy is obtained when the H atom is located close to 
the octahedral interstitial position next to an Al, Ga or 
In atom. This formation energy is compared with the 
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formation energy for two non-interacting impurities, ob- 
tained by adding the formation energy of a single intersti- 
tial H atom (Hi) and a single Al, Ga or In dopant in ZnO 
(Xzn)- The results are shown in Fig. [TJ Fig. [2] and Fig. [3] 
for the co-doping cases H-Al, H-Ga and H-In, respec- 
tively. The shown curves are for the oxygen poor limit. 
For the oxygen rich limit, all curves shift down with 3.43 
eV. This choice has no influence on our conclusions, as 
they are only based on the crossing points between the 
formation energy curves for different charge states. As 
both H and (Al, Ga or In) can give away a single elec- 
tron, the +2 charge state is the ground state for small 
values of the Fermi energy. (The +1 case is not shown, 
as it never becomes the ground state.) It is clear that 
the case of two non-interacting, thus infinitely separated, 
impurities has the lowest formation energy. In this +2 
charge state, both shallow donors repel each other. More 
importantly, however, for Fermi energy values close to 
the band gap of 3.4 eV, the neutral H-(Al,Ga,In) com- 
plex has the lowest formation energy. Furthermore, the 
crossing between the non-interacting charge case +2 and 
the neutral complex occurs in the band gap, i.e. a deep 
level is formed. This crossing occurs at 0.18 eV, 0.47 eV 
and 0.67 eV below the bottom of the conduction band 
for Al-H, Ga-H, and In-H, respectively. This is a conse- 
quence of the large decrease in formation energy of the 
neutral complex in comparison to the formation energy 
of two non-interacting neutral impurities: 1.31 eV for 
Al-H, 1.42 eV for Ga-H, and 2.11 eV for In-H. As a con- 
sequence, the shallow donor H binds to and passivates 
the shallow donor (Al, Ga or In). 

The formation of a deep level is also clear from the 
density of states (DOS). As an example, the DOS for 
the Al-H co-doped case obtained within the LDA ap- 
proximation is shown in Fig. 3J The neutral complex is 
considered in its ground state configuration. Comparing 
Figs. SJa) and (b) for the pure ZnO crystal and the co- 
doped case, a shark peak appears in the ZnO band gap. 
From the projected (local) DOS plots on H (Fig. E{c)) 
and Al (Fig. HJd)) we observe that this deep donor level 
has H s and Al p character. While an individual intersti- 
tial Hi and a substitutional Alz n atom in ZnO have both 
a defect level in the conduction band with these charac- 
ters, a bonding level is formed in the case they form a 
complex, which lowers the energy to a value in the band 
gap. 

To gain a better understanding of the physics behind 
this finding, we considered the bond length in the neu- 
tral complex. The bond between both shallow donors 
is clearly reflected in the interatomic distances. In the 
neutral case, the equilibrium distance between the inter- 
stitial H atom and the substitutional group III atom is 
1.86 A for the Al case, 1.77 A for the Ga case and 1.87 A 
for the In case. However, in the +2 case, these distances 
increase till close to 3.15 A for all three cases, showing 
a clear repulsion. As an example, for the Al-H co-doped 
ZnO, the structure of the neutral Al-H complex is shown 
in Fig. [5fa) and (b), and for the repelling +2 case in 



Fig. Efc) and (d). It is also interesting to compare these 
bond lengths in the case of the neutral complex with 
those of the corresponding molecular species A1H, GaH 
and InH, which are 1.65 A, 1.66 A, and 1.84 A 20 . The 
Al-H and Ga-H bond length are thus longer in compar- 
ison to the length in the molecular species, especially in 
the Al-H case, while the In-H bond length is much closer 
to the molecular structure InH. This shows that there 
is indeed a bond in the neutral H-X complex, although 
weaker than in the molecular case. 

To support this view and to propose an experimental 
way to detect these complexes, we also studied a vibra- 
tional property. The microscopic geometry for hydrogen 
in semiconductors is usually determined through obser- 
vation of the infrared spectrum of the related local vi- 
brational mode. In order to facilitate the experimental 
observation of the H-(A1, Ga, In) complex, we calculate 
the H-X stretching mode. As the H atom is so light we 
also take into account the anharmonic effect. Therefore 
the potential curve is fitted by the following formula^ 



V(x) = -x 2 + ax 3 + (3x A . 
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We have taken displacements up to ± 30 % of the X- 
H bond length along the bond direction. By applying 
perturbation theory to the one-dimensional Schrodinger 
equation an approximated analytical solution in the case 
of the anharmonic potential [eq. (pQ)] is given by 22 
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where uj is the frequency of the mode, ujh and Auja its 
harmonic and anharmonic contribution, and \i is the re- 
duced mass. The obtained results for the frequency of 
the local vibrational modes are shown in Table [U Again 
we can compare with the corresponding frequencies of 
the molecular species. These results are also included in 
Table U under u mo i (taken from Ref. 20 ). Note that the 
obtained frequencies for the complex are smaller than in 
their molecular equivalents. The difference is most pro- 
nounced for Al-H, and the smallest for In-H. The larger 
the difference, the weaker the bond is in the complex in 
comparison with the molecular equivalent. This trend is 
consistent with the trend that was observed for the bond 
lengths. This then also explains why the defect level of 
the complex is deepest in the In-H case. 

To summarize, we have studied by first principles cal- 
culations the co-doping of ZnO with hydrogen and a 
group III element (Al, Ga, or In). While these impu- 
rities are known to be shallow donors in ZnO, and are 
therefore expected to repel each other, we show on the 
contrary that they attract each other, resulting in a neu- 
tral complex which forms a deep level. Hydrogen thus 
passivates the group III impurity in ZnO, which puts a 
limit on its n-type conductivity. In semiconductors in 
which hydrogen shows amphoteric behavior, passivation 
of a shallow impurity does not mean that a complex is 
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FIG. 1. Formation energy as function of Fermi energy for the 
H — Alzn complex (full black curve) and for two noninteracting 
interstitial hydrogen and substitutional Alz n impurities 
(dashed red curve), in the oxygen poor limit. 



TABLE I. Calculated harmonic frequencies (cj#), anharmonic 
contributions (ooa) and total frequencies (cj) of the local vi- 
brational mode of the X-H complex, and the corresponding 
frequency {uJmoi) in the corresponding molecular structure, 
taken from Ref. 20 . All frequencies are expressed in cm -1 . 
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formed, the formation of such a complex may just be 
a consequence of the passivation. In ZnO however, in 
which hydrogen acts as a shallow donor, it is the pairing 
between hydrogen and the substitutional shallow donor 
into a neutral complex that realizes the passivation. 
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FIG. 2. As Fig.[Q but for Ga-H co-doped ZnO. 
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FIG. 3. As Fig. [B but for In-H co-doped ZnO. 
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FIG. 4. Calculated DOS, with zero energy taken at the Fermi 
level: (a) total DOS for the pure ZnO crystal; (b) total DOS 
for the Al-H co-doped ZnO, the inset is a zoom around the 
Fermi level; (c) projected (local) DOS on the H site; (d) pro- 
jected (local) DOS on Al site. 
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FIG. 5. (a) and (b) The Al-H neutral complex, (c) and (d) 
The +2 charge case of the Al-H complex. Red atom is Zn, 
grey atom is O, green atom is Al, white atom is H. These 
figures are generated by VMD version 1.8.7—. 
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